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ABSTRACT 


The intensity and direction of the geomagnetic field 
vector not only vary from place to place on the surface of 
the Earth, but they also show variation with time. The long 
term variations due to causes within the Earth are referred 
to as secular variation. The rather brief historical records 
can be extended backwards by sampling sequences of dated 
sediments, In this study paleomagnetic specimens were 
collected from a 7 metre thick section representing some 
9,000 years, which is part of a much thicker seguence of 
strata laid down during the Olympia Interglaciation (43,000- 
20,000 yrs.B.eP.)- The specimens carry a stable remanent 
magnetization, and the overall mean direction is D=005.79, 
T=+65.49 (N=94, R=93.3099, k=135, alpha95=1.3°), which is 
Significantly different from the present geomagnetic field 
vector at the sampling locality (D=022.39, I=+72.6°9). The 
mean of the virtual geomagnetic poles lies at 0.0°F and 
85.6°N (N=94, R=92.7032,K=72,A95=1.7°) and is statistically 
distinct from the Earth's spin axis being both ‘right- 
handed* and ‘far-sided'. There is no evidence of the 
recently reported Mono Lake and Lake Mungo ‘excursions’ even 
though the age of the above events are spanned by the 
sequence. Instead a sequence of smooth declination and 


inclination oscillations are observed with ma jor 
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periodicities of a few thousand years. These trace out a 
remarkable pattern in which three-quarters of an open loop 
is first traversed counterclockwise, then clockwise, and 
finally counterclockwise again. The most likely explanation 
involves an easterly biassed non-dipole field and a 
geocentric dipole which undergoes a linear feci7 Ai on with 


a period of about 5,000 years. 
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CHAPTER 1 


INTRODUCTION 
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1.1 Geomagnetism and Secular Variation 


The Chinese, probably as early as the second century 
BeCs, invented the earliest known form of magnetic compass 
(Needhan, 1962). In Europe the directional properties of the 
lodestone were referred to by Alexander Neckham in 1190, and 
Petrus Peregrinus discovered the dipole nature of the magnet 
in 1269. The magnetic declination was known to the Chinese 
in the eighth century A.D. (Needham, 1962; Smith & Needhan, 
1967), and was rediscovered in Europe in the latter part of 
the fifteenth century. The magnetic inclination was 
discovered independently by Georg Hartman in 1544 and Robert 
Norman in 1576. In 1600, William Gilbert synthesized these 
observations in his well known treatise ‘De Magnete’,. 
Gilbert described the Earth's magnetic field as that of a 
uniformly magnetized sphere. By 1839 global coverage was 
sufficient for C.F.Gauss to undertake the first thorough 
Mathematical analysis of the field. His spherical harmonic 
analysis showed that more than 99.5% of the Earth's magnetic 
field is internal in origin and all subsequent analyses have 
confirmed this finding. This type of analysis further shows 


that approximately 80% “of ‘the Earth's field (can be 
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attributed to a single geocentric dipole inclined at 11.5° 
to the axis of rotation with a magnetic moment of about 
8x1022 As.m2. Approximately 85% of this axial dipole field is 
attributed to a single dipole along the rotation axis and 
the remaining 15% to two orthogonal dipoles lying in the 
equatorial plane. When the best fitting dipole field is 
Subtracted from that observed over the Earth's surface the 
residual is termed the non-dipole field. At present there 
are some eight non-dipole features (foci or anomalies) of 
continental dimensions displaying either positive or 
negative values with typical amplitudes of 10-S T (Fig. 
1.21). Lowes (1955) and Alldredge & Hurwitz (1964) 
represented the dipole and non-dipole fields by a central 
dipole and an array of 8 to 10 smaller eccentric radial 
dipoles with moments ranging from less than 1x1022 A.m? up 


to 5x1022 A.m.2. 


As yet there is no completely satisfactory answer to 
the problem of the origin of the Earth's magnetic field. The 
most plausible theories involve some form of dynamo action 
(Elasser, 1946; Bullard, 1949). The Earth's core is both a 
fluid and a good electrical conductor. In other words it 
permits both the flow of electrical current and mechanical 
motions: the resulting interactions are thought to generate 
a self sustaining magnetic field. A comprehensive recent 


review is given by Gubbins (1974). 
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Fig. 1.1 The non-dipole field for 1945. The arrows indicate the 
direction and magnitude of the horizontal component. 


The contours give the vertical component at 


intervals of 2.0x107°T (after Bullard et al., 91959); 
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The intensity and direction of the magnetic field 
vector not only vary from place to place on the Earth's 
surface but they also show variation with time. There are 
two distinct types of temporal change; transient 
fluctuations and long term variations. The transient 
fluctuations arise from causes outside the Earth; aes are 
very small and produce no enduring changes in the 
qeomatneeie field. They are of no consequence in the present 
study. On the other hand the long term variations, often 
referred to rather loosely as secular variation, are of 
prime interest to us. These are due to causes within the 
Earth and the net effect over long intervals of time may be 
considerable. The ‘spectrum*® of temporal variations of 
internal origin exhibits periods ranging from less than 102 


years to 106 years or more (Table 1.1). 


At the surface of the Earth the non-dipole component of 
the field, although much weaker than the dipole, changes 
more rapidly, with periods as short as decades. There seems 
to have been a general westward drift of approximately 0.29 
of longitude per year, over the first half of this century 
(Bullard et ale, 1950) but paleomagnetic studies have 
indicated the presence of both eastward as well as westward 
drifts (Yukutake, 1962; Skiles, 1970). Examples of 


historically recorded changes in declination and inclination 
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Table 1.1 The Spectrum of Temporal Variations of the Geomagnetic 
Field of Internal Origin ( after Jacobs, 1975 


McElhinny & Merill, 1975 ). 
Types of variation 
Minimum observable period 


Secular variation of non-dipole field 
"High frequency’ oscillations 
"Medium frequency’ oscillations 
"Westward drift' 


Secular variation of dipole field 
Dipole ‘wobble’ (eastward drift 
Field strength oscillations 


Main field generation mechanism 
Magnetic Excursions 


Polarity Reversals 
Change in field direction 
Change in intensity 
Time interval between polarity changes 
Bias in favour of one polarity 


Period 
B87 
<1100 


100 — 5,000 
e000 


15200 — 1,800 
9,000 

10342 "10" 

103 — 104 
1,000 — 4,000 
3,000 — 10,000 


102 = 107 
108 = 10? 


and 


years 


years 
years 
years 


years 


years 


years 


years 


years 
years 
years 
years 
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at several observatories are illustrated in Fig. 1.2. These 
patterns usually take the form of open loops, and Runcorn 
(1959) and Skiles (1970) have shown that clockwise looping 
results from westward drift of non-dipole features and vice 


verlTsae 


One of the most Significant discoveries in 
paleomagnetic studies is that of reversals of the Earth's 
magnetic field. In the dynamo theory a field could be 
produced in either direction. There is no a priori reason 
why the magnetic field of the Earth should have a particular 
polarity and there is no fundamental reason why its polarity 
should not change. The observed field reversals are highly 
irregular and therefore have been modelled aS ae non- 
stationary random process. Cox (1968) developed a 
probabilistic model to explain reversals in which it is 
assumed that the changes in polarity result from an 
interaction between random processes and steady 
oscillations. The steady oscillator is the dipole component 
(Fig. 1.3) whilst the random variations are due to the non- 
dipole field. The random variations serve as the triggering 
mechanism that results in a reversal whenever the non-dipole 
to dipole ratio exceeds a critical value (Fig. 1.4). Cox 
(1975) suggests that the non-dipole sources are not randomly 
distributed with respect to latitude. During times of normal 


polarity, flux emerges from the core in localized regions 
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Fig. 1.2 Secular change of declination and inclination 
at 


Baltimore, Boston, London and Paris from 
observatory measurements. Note the same 
clockwise motions even though the size 

the loops vary (after Nelson et 
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DIPOLE moment (1077 A -m” ) 


sd 
) 1000 2000 3000 4000 5000 6000 7000 8000 9000 
AGE (YEARS BEFORE PRESENT) 


Fic. 13 gVariatdon. odin geomagnetic dipole moment from ° 
archaeomagnetic studies. The short bar near the 
ordinate represents change during the past 130 
years as determined from spherical harmonic 
analysis. The archaeomagnetic results are averages 
over 500 year intervals. The number of averaged 
data is shown above each point and the vertical 
lines indicate the standard error oof the mean. 
Unshaded squares refer to intervals containing 
too few data to provide meaningful statistics 


( after Cox, 1968). 
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Fig. 1.4 Probabilistic’ model for reversals. T is the 
the dipole field and- 7” ‘the dlensth of a 
interval. Whenever the quantity (Ma + Ma" ) 
Bien <a. POLArLCy change - witt occur. Ma, is 


moment of the dipole field and Ma' is the 


non-dipole moment (“after ~Coxtie 1900"). 
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near the equator and reenters the core in localized regions 
approximately 55° from either pole. On the other hand during 
times of reversed polarity this entire pattern is reversed 
suggesting that the toroidal field also reverses when the 
dipole field reverses itself. The latitudes at which the 
large features of the non-dipole field are generated by 
convective cells or fluid motions therefore appear to be 


linked to the polarity of the dynamo. 


Apart from reversals, the Earth's field appears to 
undergo ‘excursions’ which seem to represent a kind of 
enhanced secular variation but are not total reversals. 
Barbetti & McElhinny (1976) give a summary of localities and 
ages of geomagnetic excursions recorded over the past 40,000 
years. However in a recent review Verosub & Banerjee 
(1977,p.153) “do not feel that the existence of any proposed 
excursion is yet sufficiently well established for its 


reality to be beyond question." 


1.2 Archacomagnetism_and Paleomagnetisn 


The records of long term variations (>100 years) of the 
Earth's field can be extended backwards for a few thousand 
years using remanence preserved in historically dated lavas 
and archaeological material. The main difficulty encountered 


is the precise dating of the samples, the intermittent 
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extrusion of lava flows and the sporadic rise and fall of 
past civilizations. An alternative method is to sample 
sequences of dated sediments, especially in the Late 
Quaternary where the C1* method of age determinations can be 
used. The geomagnetic record in these sediments is primarily 
acquired by a process known as detrital remanent 
magnetization (DRM) (for a recent review see Verosub, 1977). 
Renee tonal DRM is the process whereby magnetic particles 
eroded from preexisting rocks are aligned by the ambient 
field as they settle. The process of rotation of such 
particles into the ambient field direction when they are in 
the water-filled interstitial holes of a wet sediment prior 
to consolidation is termed post-depositional DRM. Johnson et 
ale (1948), King (1955) and Griffiths et al. (1960) have 
investigated the depositional process under controlled 
laboratory conditions whilst Irving & Major (1964) and Blow 
& Hamilton (1977) have investigated the post-depositional 
process. The laboratory studies of post-depositional DRM 
show impressive agreement between the direction of the 
applied field and that of the magnetization of the sediment, 
however studies of depositional DRM often indicate an 
‘inclination error’, iee. even though the declination is 
faithfully recorded by the sediment, the sediment 
inclination is less Oona n the applied field inclination. 


Irving (1957, 1967) and Opdyke (1961) have confirmed that 
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these effects do not seem significant for most rocks of 
paleomagnetic interest since any inclination error acquired 
during the deposition of the particles is ‘corrected for’ by 


Subsequent post-depositional DRM. 


1.3 Organization of Thesis 


The present work is a detailed follow up study of a 
reconnaisance investigation by C.J. Oberg (1978) . A 
Quaternary section which spans some 9,000 years was sampled 
at stratigraphic intervals representing approximately 100 
years to achieve a high resolution data base suitable for 
investigation of secular variation patterns and application 


of spectral analysis techniques. 
The subject matter is divided as follows: 


Chapter 2 gives an account of the field and laboratory 
techniques employed. The rock magnetic froperties and 
mineralogy of the magnetic carriers are discussed in Chapter 
3. Analysis of the paleomagnetic results is discussed in 
Chapter 4. The maximum entropy method is used to obtain 
improved knowledge of the geomagnetic spectrum since it is a 
radically new approach to power spectrum estimates which 
obviates several limitations of conventional methods and is 
especially suitable for short time series. Worldwide 


geomagnetic secular variation results during the Late 
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Quaternary are summarized in relation to the present study 


in Chapter 5. 


fhe standard Statistical techniques due to Fisher 
(1953) are used in analysing the data presented. Fisherian 
Statistics are essentially an analogue of Gaussian 
Statistics adapted to distributions on a sphere. The most 
important parameter involved is the precision parameter, 
k= (N-1)/(N-R) ( N = number of unit vectors and R = length of 
the resultant of the N unit vectors ) which corresponds to 
the Gaussian invariance. At a probability level of (1-P), 
the true mean direction of the population lies within a 
circular cone about the vector resultant R with a semi-angle 


1/(N=1) 


alpha =cos-1| + SF { (1/2) aoe Hop ssmaliw alphapwithe 


C1—P) 
circle of 95% confidence is given by alpha), = 140/(KN) 172, 


Se-Ie units are used throughout the thesis. Since much 
geophysical literature is still written in c¢c.g.S units, 


sonversions for the most common parameters are given below. 
Magnetic dipole moment, 1 A.m@ = 103 emu 
Intensity (moment/unit mass), 1 A.m@/kg = 1 emu/gnm 


Magnetic induction, 10-9 Tesla = 10-5 gauss = 1 gamma 
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CHAPTER 2 


Detailed Paleomagnetic Results froma 


Quaternary Section in British Columbia 


' The late Quaternary geologic/climatic units in southern 
British Columbia and northwestern Washington State record a 
sequence of glacial and interglacial episodes, usually 
tezerred to as Postglacial (10,000-0 yrs.B.P.), Fraser 
Glaciation (20,000-10,000 yrs.B.P.), Olympia Interglaciation 
(43,000-20,000 yrs.B.eP.) and Okanagan Glaciation (>43,000 
¥rs.B.P.) (Westgate and Fulton, 29S) = The Olympia 
Intergiacial sediments in the interior plateau region of 
British Columbia contain several thin layers of rhyolitic 
and dacitic tephras which represent the distal portions of 
widespread airfall eruptive units. Each tephra layer, being 
of regional extent and distinctive Al character, 


serves as a useful time-parallel stratigraphic marker. 


In this study paleomagnetic specimens were coliected 
from a 7 metre thick section which itself is part of a much 
thicker sequence of strata laid down during the Olympia 
Interglaciation (Fig. 2.1). The outcrop sampled forms a 


stream-cut escarpment on the edge of Bessette Creek 
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Stratigraphic succession at the Riggins Road 
exposure. The section sampled is part OF 
a sequence laid down during the Olympia 
Interglaciation 43 000m— OR OO00 vrs. Babe). 
The four radiocarbon dates are shown. 


15 


ae 


i y pene 
ie , tet v 
ntebe Rie 
Tf Fed) 2; 
Bea ea 
Se \ ah 
Mea sr 7 : 


“a 


16 


Maes r (330¢q ) 
dg °''006FOOZIE 1€07-DSO 


w0upag 


S}IBOGED ujo}d500}} Sozudo, puo ipad 
*Ppoom ‘sxsnow usyim fyIS puo "puds Ysissse 620 


¥[ ts 9: 
PUBS D/0}}09 pud {Is duUlsysNdD]-0/9019 
“puos Aigqoa Ne36 hia 
5 = “p pus jaa0s6 Kiqqo5 
02d 9 55408 7 F OOLLe scél Se Bk IS puo puos 
oSeeisINta cs! gDad [Oylssap fy yyy oC, *}8a045 50 $asud) jod0) Yyim {tL 
4yis 6 WHS @UIGSNIO{-0190)9 


TI Dsydsay POOY suibGhiyiststs QN3937 


4 
¢ 
ahs © 
E 
ec 9 
: ( poom aodnids ) S 
_- 4 8 S407TEFOOESZ ES61=9S9 a 
ete < saho] yous DiunB10' 4, Ne 7). © 


i 
4aApj 4418 Apups aut 


apts 


l 
paq Apups BSIDO IST: 
4]! ie” es 
, (309d 41449p) _ 
re hous Syél-> 
en TSBHMOngROOTST, EAI Oe S. 


F 


vere SONDINO GB sapuy 
Py OIC WAIO 


UO}}0190/9 
408044 


<0) 


nn ARES. 


17 


approximately 10 km east of Lumby, BoC. (lat. 50°7'55"N, 
long. 118951°45"W), and consists of dark organic-rich flood 
plain sediments. With the exception of a coarse sandy bed 
approximately 0.3 m thick, and a fine sandy silt layer 
approximately 0.2 m thick near the top of the section, the 
sediment is remarkably uniform in character and there is no 
evidence of any major change in depositional environment 
during the time span sampled. The section contains four 
radiocarbon dated horizons and two tephras, namely the 
Cherryville and Riggins Road tephras (Westgate and Fulton, 
1975). The radicmetric dates are vital for estimating the 
average rate of deposition, which assumes considerable 
importance in later analysis and discussion. Two of the 
dates were obtained from wood specimens, but the other two 
were from detrital peat. Some time may have passed between 
the death and deposition of any biological material, 
therefore both wood and detritus dates may be too old for 
their stratigraphic position, but they cannot be too young. 
The wood, having retained its integrity, is less likely to 
have undergone a long pre-depositional history. Since both 
detritus dates appear too old relative to the wood dates 
(the upper one is even out of sequence), the sedimentation 
rate given by the two wood dates is likely to be the most 
accurate. The lower wood date (31,200+900 yrs.B.eP.) lies at 


the base of the section. The upper wood date (25, 3004320 
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yoSseBeP.) lies 4.50 m above this, giving an average 
deposition ‘slowness' of 1342 yrs/cm (deposition rate of 
0.77+40.01 mm/yr). The discrepancies between the various 
dates are not large - a linear regression using all four 
dates yields a deposition ‘slowness! of 10 yrs/fcm. The age 
of the oldest specimens is nominally 31,200 yrs.B.P., and 
extrapolation of the estimated deposition rate yields an age 
of 22,100 yrs.B.Ps for the youngest specimens coliected. 
Thus the total time span represented is slightly over 9,000 


years. 


2.2 Field _and Laboratory Techniques 


The specimens were collected in commercially available 
plastic cubes of 2 cm edge. Two small holes were drilled 
close to the base of each cube to let out air whiist the 
cube was being pressed or tapped gently into the vertical 
face of the outcrop. Each cube was painted on one face in a 
way which permitted unique orientation. Specimens were 
oriented by measuring three angles, strike, dip and roll 
with a Brunton compass-clinometer, as described in detail by 
Oberg (1978). Over the course of six days sampling, forty 
magnetic sun bearings were taken giving an average local 
magnetic declination of 21.7°E+40.8°9 (s.sd.). These data allow 
the orientation scheme to be assessed. The mean agrees guite 


well with the regional average (22.3°E, Canadian Isogonic 
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Chart for 1975.0), indicating that adequate accuracy was 
achieved. Likewise, the smail variance, of less than a 


degree, indicates that good precision was attained. 


After removal from the outcrop, each specimen was 
capped and sealed to prevent dessication. A total of ninety 
four sites (i.e. stratigraphic horizons) with four 
independently oriented specimens per site was sampled, 
beginning immediately below the Cherryville tephra. Sample 
Spacing was as uniform as the outcrop allowed; of the ninety 
three intervals between sites, eighty nine were 7 cm, two 
were 9 cm, one was 12 cm and one was 18 cm. The 
sedimentation ‘slowness’  (134+2 yrs/cm) indicates that the 
sampling interval is about 100 years, with each specimen 
representing 20-30 years. At each site the lateral spread 
was typically 30 cm but ranged from 10 to 60 cm. All 
Sepa aeertae elevations in this thesis are measured 
downwards from the Cherryville tephra (Fig. 2.1). The layer 
which extends from 0.63 m to 0.93 m was sandy and tended to 
crumble upon sampling, however the occurrence of a thin silt 
layer within this sandy layer allowed us to avoid too large 
a sampling gap. A total of eight sites were sampled within 
the hard compact silt layer which extends from 4.40 to 4.89 
Me These specimens were obtained by carving stumps 
protruding from the vertical face of the outcrop over which 


plastic cubes could be snugly fitted. A comparison of data 


et + | haa 


au 42 SIOBOS . ate pops tan odineltad * (Osave r aon d 
& tafe feasl la ‘ Vennerisy fisse ) ont yaedwo4t /- 
Pan ft6S7R TEV Ao Fe LOossG: DOOD iia saxsoioas 


Be 
ES 


cw nebinece dose doo THO od t AGae Tpsguet aaa 3A 


Pe 


Vioain Bo [edn 4 «hOLtesiceaat toavenq oF hetsee ‘bate he iq 
tbe Av ty janes tion > mm sipiteate sith  savse 
v—betease 254 , oti g en *m£2Sqe Paaerpee . ¥E sie oar: 
yigese. .ai1tge? oi ‘ va | vat ‘woist viststiguns osiae 


vaegin of) i¢} shawoll!s Tipichi ls ait 2h ‘gtotinan as sl pakoa 


awd 7 g 7 ; lav & tia ¢tuini Aes tLe, usevted es ¥oaetvar 4 
( mi i 
ivy ef Ad 26 } a0 i 5 ie Ta cr 25K, : ane ; e8D 
aie ved? - ray Leg Sees. Fase aioia? OF FHS 
ia a ; ‘AE 
aeuioege: isso Wiiw yemsdy OG: Sods =f bees eaai gailgue 


Pte 
> 


ait 39.te Haas shal Set ON F QE=0; Spada’ 

sao- OW | Ot FOF «0072 | bapa x Jud -md OE yLisog 

uahs se ‘5  akesed? bet at wanroitayele ise 
ayek. “eat iff 5) ote) tea Ca ae bdt w072 4 
of ‘bebast bn vhage 25W m Ee a0" OR a ea.0 O43 afinst x9, 


nfd? 4 Bersonea3es So + ots spans spe hga es aoqa 


+ a 


Marcel "ee brows oF Hp po venite ‘age, tems Bis nines rome 
dissis Dalla ese STB a3 he tips te esos & +450 et aR 
BULO- - oF 00.1) Ox? 2h06¢%9 dokde ways ahs AAAS bish one 


or 


ayo? Le & bee oo) “od He 0.0534 6 osu * anateSoge seog? 7 
foivs moro’ qoassuo edt bc epee Psbit2sv att ieee chain oy) 


"i ; 


Blah FO Ace pgn00 A alata teudéae bfuoo “zeta sbrbeig 
ios) 


— : | * ao 


20 


from *carved* versus ‘'tapped' specimens indicates that any 
distortion taking place as the cube penetrates the sediment 
is minimal and its magnetic effects negligible. The data 
involved in this test are brought together and discussed in 


section 2.4, 


The specimens were stored in a magnetically shielded 
room (field strength approximately —40nT, Reid, 1972) and 
taken out only for remanence measurement using a Digico 
balanced fluxgate magnetometer. All measurements were 
carried out with an integrating time of 37 seconds (i.e. 284 
spins), and instrument calibration was repeated every six 
measurements (ie¢€. approximately every 30 minutes). Empty 
cubes, painted and unpainted, gave no detectable signal 
above the noise level expected when the spinner is operated 
with no sample present (Molyneux, 1971). The specimens 
themselves are at least two, and generally three, orders of 


magnitude stronger than this. 


Schmidt (1976) constructed Helmholtz coils to cancel 
the effect of flux concentration at the top of the mumetal 
shield of the magnetometer, which amounts to about twice the 
Earth's field. This potential problem was investigated by 
exposing a typical specimen to a field twice that of the 
Earth (approximately 0.1 mf) for periods from 10 minutes to 


4 hours. No detectable remagnetization was observed and 
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therefare the flux concentration claimed by Schmidt does not 
seem to be a problem in this study. Perhaps Schmidt suffered 


from strong but very soft and unstable material. 


Routine stepwise alternating field demagnetization 
studies were undertaken using the apparatus described by 
Murthy (1969) which is essentially identical to that 


developed by McElhinny (1966). 


2.3 Results 


To attempt isolation of the primary magnetization from 
magnetic *noise', incremental demagnetization was 
undertaken. It is both impracticable and unnecessary to 
carry out complete demagnetization analysis of every 
specimen in‘an extensive collection such as this, instead 
the routine procedure of investigating selected pilot 
Specimens in detail was followed. The treatment which best 
isolated the stable directions was then applied to the 
remainder of the collection. Ten pilot specimens were 
selected (one each from sites HAA 0 m, HAK 0.75 m, HAO 1.07 
m, HBF 2.28 MN» HBI 2.49 m, HBT 32.26ime HCE19.05¢cm,0HCO, 4Eg5 
m, HDD 5.80 m and HDL 6.36 m) and subjected to stepwise AF 
demagnetization in peak fields of 2.5, 5.0, 10, 15, 20, 30, 
40, 50 and 60 mT At this point three specimens (HAK02, 


HAOO3, HBTO4) remained stable and were therefore treated in 
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higher fields; however they all underwent large angular 
Shifts at 80 mI. Median destructive fields (MDF's) range 
from 19 to 36 mT with an average close to 20 mT (Fig. 2.2). 
Up to 30 mi, successive directional shifts were small 
(mean=2°, maximum=7°) ; beyond 30 mT some, but not all, 
specimens became unstable (Fig. 2.3). The results are 
illustrated by the Zijderveld (1964) convention in Fig. 
2e4a-j. Also plotted is the rate of angular shift between 
successive demagnetization steps which has been used by some 
authors as a measure of stability (Symons & Stupavsky, 1974; 
Lowrie & Alvarez, 1977). As a result of the pilot 
investigation 10 and 20 mf were chosen as appropriate fields 
for AF cleaning of the remaining 366 specimens. Again, 
angular movements are small (mean=2° , for both the 0 to 10 
mT and 10 to 20 mT steps), and the overall means at 10 af 
and 20 mf are not significantly different. It thus appears 
that these samples carry a stable ‘remanence, whose origin 
and significance is the subject of the following chapters. 
Prior to cleaning, intensities ranged from 2.94x10-6& to 
4.28x10-4 Aem2/kg with a geometric mean of 2.77x10-S 
A.em2/kg. After cleaning in peak fields of 10 and 20 mT, the 
geometric mean intensities are 2.12x10-5 and 1. 45x10-5 
A. m2/kg he sset wells. An archive of all the data (1,208 
individual determinations of remanence) is given in Appendix 


1. The results are illustrated in terms of magnetograms in 
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Fige 2e5a-c, which show all the individual specimen data at 
0, 10 and 20 mf. It is clear that smooth declination and 
inclination oscillations are present and the significance of 


these will be discussed in chapters 4 to 5. 


2-4 Comparison of Sampling Methods 


As discussed above, two sampling techniques were 
employed. The majority of specimens were collected by 
tapping the plastic cubes directly into the face of the 
Outcrop, but where the sediments proved too hard carving was 
used. We also collected a large coherent block spanning the 
interval 1.28 m to 1.40 m. The results obtained by these 


different techniques are now assessed. 


From the large block a total of five sites (four 
specimens each) were carved in the laboratory. The 
declination and inclination of these five sites in relation 
to adjacent collecting horizons are shown in Fig. 2.6. It is 
clear that the two data sets are in clcese agreement. 
Orientation error of the block, whose frontal face is not 
perfectly planar, is the most probable cause of the slight 
systematic differences observed. The elevation of one of 
these five sites ( HEG ) corresponds exactly to site HAS. 
Two other horizons ( HEB = HBE, elevation 2.21 m; HEA = HDC, 


elevation 5.73 m ) were sampled by both methods in the 
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Fig. 2.6 Site mean results from laboratory carved sites 


and 
the corresponding sites sampled in situ. Tapped 


Specimens were demagnetized at 10 mT. The large 
coherent block from which the carved specimens 
were obtained was not stored in the magnetically 
shielded room. The end points of these specimens 
were established at 20 nT, except for specimens 
from site HEG where the end points were at 35 nT 
( see™ footnote to Tablé 2.1). "The “bars indicate 
6D and 6I where 6D = alphag./cosI and 61 = alphags . 
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field. It is therefore instructive to quantitatively compare 
results from these three elevations. Table 2.1 gives the 
appropriate data; again it is clear that differences are 
very small. The standard F-ratio test (Watson & irving, 
1957) indicates that these small differences are not 
Statistically significant. The data of all the carved 
Specimens are given in Appendix 2. One further comparison 
between the sampling techniques is provided by the within- 
Site dispersions. For the 86 ‘tapped* sites the precision 
parameter (k) ranges from 20 to 5,242 with a mean of 342 
Whilst for the 8 ‘carved* sites it ranges from 140 to 568 


with a mean of 274 (all data at 10 mT). 


A point of considerable geomagnetic significance 
emerges from the data illustrated in Fig. 2.6. The finer 
Sampling scale reveals no prominent geomagnetic features and 
this implies that the 7 cm spacing used is generally 
adequate and that no persistent high frequency content is 


missed. 


. 205 -Copparison .of ithe Present Stud with _Oberg‘s 


Reconnaissance Results 


A comparison is now made between the present results 
and those reperted by Oberg (1978). The distance between the 


Cherryville tephra and the top of the gravel layer is 
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Table 2.1 Comparison of sampling techniques 


Elev (m) 


= 


Carved (in field) 
EE 


HAS 


1.35 


All data refer to AF treatment at 10 mT except site HEG 


(data marked by asterisk ). The specimens from this site 
were unintentionally placed close to a strong magnet EOL, 
a period of almost one week , and thus required further 
AF treatment ..The end points were established after 


treatment “at.-35 mi =. 


A angular difference 


DONE BUR, Ry 
B- 1 DUN, -> Ry 


where B number of sites 
Ny number of specimens at éach site (i.e. 4) 
Ry; length of the resultant vector at each site 
and Rg length of the vector sum of the resultant 
at- each, site 
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reported by Oberg to be 6.91 m, whilst in this study it is 
6.71 Me Since the previous investigation was a 
reconnaissance study only (all the sampling was done in a 
Single day), the present elevation measurements are 
considerably more reliable.t A comparison of the two 
intensity profiles shows remarkable similarities, reflecting 
the same underlying lithological variations (Fig. 2.7). When 
the former data are compressed by a factor of 0.97 ({iee. 
6.71/6.91) the agreement becomes extremely good (Fig. 2.8), 
giving one much confidence in the results. The corresponding 
declination and inclination magnetograms are also shown. The 
present data exhibit oscillations of somewhat larger 
amplitude at certain elevation intervals, but on the whole, 


the agreement is satisfactory. 


2 oe @@ @ ee 2 2 @@ @ae@ a= = 


10berg and his field party agree that this is the case. 
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Fig. 2.7 


Comarison of intensity yeprofiles 
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Fig. 2.8 Comparison 


site mean intensity, 
declination 


and inclination 
study with Oberg's reconnaissance 
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CHAPTER 3 


The Nature of the Remanent Magnetization 
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In Beneoa anes te studies it is important to enquire if 
the remanence finally isolated by partial demagnetization is 
‘primary’ (ieee aligned parallel to the ambient field at the 
fine (of ormation)...The direction and magnitude of the 
remanence removed by demagnetization can yield useful 
information about the nature of any secondary magnetization. 
Simple vector subtraction yields the direction and intensity 
of the remanence removed at each step by partial 
demagnetization. Such calculations have been performed on 
the 94 sites Of. the Riggins, ykoad: data wins, the fAF 


demagnetization ranges of 0 to 10 mf and 10 to 20 mT The 


results are illustrated in Fig. 3.1a-c. 


Consider first the means of all 94 site means after 0 , 
10 and 20 mf treatment (Fig. 3.1a); these represent the 
average vectors remaining after these treatments {see 
section 4.2 for detailed discussion), not the vectors 
removed by the treatments. The 10 and 20 mf data are not 


Significantly different, but the NRM data are very slightly 


displaced towards the present field vector. This suggests 
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Fig. 3.1(a) 


(b) 


(c) 


The 


Mean of the 94 site means at 0 ; 
10 and 20 mT. (Azimuthal equidistant plot) 


Mean of the 94 site means removed 
vectors at 0 —10 mI and LO. = 20 it ae 
(Azimuthal equidistant plot) 


Removed vectors for the 10 pilot specimens. 
The triangles represent the individual 
O - 5 mI removed vectors ; three of these 
(shown as inverted triangles) are rather 
divergent . (Equal area plot) 


data plotted are as _ follows : 
D? 
Present Field 022. 
Site means O mT 006. 
10 mT 005. 
20 mT 005. 
Removed vectors 0 -— 10 mT 009. 
10 = 20 mT 004. 


Pilot specimens 


removed vectors 0 -— 5 mT (10 specimens) 018. 
5 -—10 mT (10 specimens) 004. 
O —- 5 mT (7 specimens) 021. 
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that a small, neeant Uypsacstaneae secondary magnetization is 
present - but the data by no means establish this fact. A 
clearer picture emerges when the removed vectors are 
considered: the 10-20 mf data lie close to the final cleaned 
direction, but the 0-10 mf data are clearly displaced 
towards the present field vector (Fig. 3. 1b). This feature 
also manifests itself as a small kink in the Zijderveld 
diagrans (Fig. 2.4a-j)- It appears that a secondary 
component aligned along the present field is being removed, 
but that the first demagnetization step (0 to 10 mT) also 
removes some of the primary remanence. Further scrutiny must 
rest on the 10 pilot specimens since these provide the only 
results below 10 mf. The 10 removed vectors at 5-10 mT agree 
with the previously established "primary’ result but the 0-5 
mf data lie much closer to the present field. Three of these 
removed vectors are quite divergent (Fig. 3.1c), but a more 
tightly grouped subset of 7 vectors is in excellent 
agreement with the present field vector. Unfortunately with 
only 10 pilot specimens involved the circles of 95% 
confidence are rather large, but there can be littie doubt 
that a recently acquired remanence is present in this 
material and that peak alternating fields between 5 and 10 
mT completely erase this and reveal an earlier, more Stable, 


and presumably primary remanence. 
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As in many other geophysical and geological problems, 
time is an important factor. Exposure of magnetic material 
to magnetic fields over long periods of time can produce 
serious effects, in particular the gradual acquisition of 
viscous remanent magnetization (VRM). It has been shown that 
Many Magnetic materials exposed to a magnetic field will 
slowly acquire a magnetization in the direction of that 
field. The process is temperature-dependent and also depends 
logarithmically on time. This is attributed to -randon 
thermal fluctuations causing irreversible magnetic changes 
and thus leading to a greater magnetization in the direction 
of the applied field. The process of VRM is important in the 
present investigation because it was found that similar 
specimens from the reconnaissance study (Oberg, 1978) 
Changed their magnetization by small, but detectable, 
amounts when stored in the laboratory over a period of two 
years. Consequently it was decided to quantitatively 
investigate the VRM properties of this material. Four of the 
pilot specimens (HAA03, HAKO2, HBIO1 and HCOO4) were first 
demagnetized in a peak alternating field of 180 mT to 
destroy any remaining NRM. The specimens were then placed in 
the magnetically shielded room and a field of 1 mT was 
applied in a known direction. The VRM's were monitored at 


suitable intervals for 100 hours. The remanent magnetization 
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of each specimen was observed to increase slowly and its 
direction approached that of the applied field (Fig. 3.2).. 
The VRM acquired initially interacted with the small 
magnetic vector remaining after 180 mT treatment but after 4 
cre the growing VRM became dominant, and the magnetization 
(J) increased logarithmically with time (t), which is 
Characteristic of VRM (Shimizu, 1960). Thus, J=S.log(t), 
where S is called the magnetic viscosity coefficient. The 
appropriate value for each specimen was determined from the 
Slope of the corresponding regression line from 4 to 100 
hours. These values of S can be used to obtain the value of 
the VRM that could have been acquired by the sediments in 
the Earth's ambient field of approximately 0.06 mT (field at 
lat. 509N) in about 30,000 years (Table 3.1). These VRM's 
are not negligible compared to the original NRM. However, 
the viscous remanence, as expected, is much softer than the 
NRM, with MDF's between 2.5 and 10 mT compared to values 
around 20 mT for the NRM (c.f. Figs. 2.2 and 3.3). This 
means that as AF treatment proceeds the influence of any VRM 
is continually reduced (see column 6 of table 3.1). The VRM 


data is given in Appendix 3. 


It is of interest to quantitatively compare the 


the magnitude of the secondary magnetic overprint removed 


from the NRM by AF treatement (section 3.1). Table 3.2 gives 
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Table 3.1 Viscous remanent magnetism of four specimens 
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@ correlation coefficient obtained from least-squares 
regression line between magnetization and 1og;, (hours) 
(for the interval 4 to 100 hours). 


b Magnetic viscosity coefficient (x 10°78 A.m? Nes 
Values quoted are those appropriate for an 
ambient field of 0.06 mT (field at lat. 50°N). 


CC ovRM (x 1078 nee ) which would be acquired 
by exposure to a field of 0.06 mT for 
30,000 years. 
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Fig. 3.3 Mean ( and range ) of the AF demagnetization 
curves of the 100 hours VRM . 
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Table 3.2 VRM's and vectors removed by AF demagnetization 


Removed Vectors Total 


specimen NRM VRM 
O10 mi ato me VRM 


All moments are quoted in A.m* x 1078 ; 
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ae appropriate data. Considering the enormous extrapolation 
involved, from 100 hours in the laboratory to 30,000 years 
in situ, the agreement is quite encouraging and there is no 
doubt that the measured VRM characteristics are capable of 
explaining the observed secondary magnetic components. 


3.3 Magnetic Mineralogy 


Any interpretation of the remanent magnetism in a suite 
of rocks depends on a knowledge of the minerals carrying the 
remanence. One test to distinguish the primary magnetic 
mineral is the study of isothermal remanent magnetization 
(IRM). The IRM of hematite and magnetite of various grain 
Sizes have been studied by Rogquet (1954). It was found that 
hematite and magnetite reach saturation at widely differing 
magnetic fields; hematite in a field of 3 T and magnetite in 
a field of 0.3 T. Four specimens (HAHO4, HAROY, HASO4Y and 
HBBO2) were first demagnetized in a peak alternating field 
of 180 mT, and then introduced between the poles of an 
electromagnet and acquired IRM in successively higher fields 
from 0 to approximately 1 T. The IRM*s were measured after 
each application of field (Appendix 4) and the resuits are 
iilustrated in Fig. 3.-4a-d. It is observed that the 
specimens saturate in a field of approximately 0.3 T. The 
IRM test clearly indicates that the magnetic carrier is 


Magnetite, (Roguet, 1954), or a related ferrite (see e.g. 
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Day, 1977). 


A second method of identification is to determine the 
Curie temperature of the magnetic constituents. The 
variation of intensity with temperature of four specimens 
(HAHO2, HARO3, HEGO5 and HBBO3) was determined with a high 
temperature Digico spinner magnetometer. The intensity at a 
particular stabilized temperature was measured over an 
intergrating time of 37 seconds (ise. 28 spins) with the 
furnace switched off. The high temperature data is given in 
Appendix 5. The behaviour of the specimens as the 
temperature was cycled from room temperature to 600°C and 
back to room temperature is shown in Fig. 3.5a-d. All four 
blocking temperature spectra are very Similar, and imply a 
Curie temperature of approximately 450°C, which corresponds 
to a titanomagnetite of composition 80% magnetite and 20% 


ulvospinel. 


It would appear that the remanence in these samples is 
carried by titanomagnetite grains, hematite being either 
absent or negligible. This, in turn, suggests that the 
remanence is depositional and is free of chemical remanence 


(CRM) overprints which generally reside in hematite. 
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Fig. 3.5a-d 


Variation of intensity with temperature. 
The intensity J is normalized to its 
value , Jpr , at room temperature . Specimen 
HEGO5 is carved aie! the laboratory from 
the large coherent block (see section 2.4). 
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CHAPTER 4 


Analysis of the Riggins Road Paleomagnetic Data 


4.1 Introduction 


In connection with the basic data illustrated in Fig. 
DeaoeovTcpe Hhebe*teare “two. topics! wf major importance to be 
discussed. These are the overall average direction of 
remanence (and its corresponding paleomagnetic pole) and the 
smooth declination/inclination oscillations. Both of these 
are intimately related to the nature of the geomagnetic 
field- and "ats’» spectrum of temporal varaations but pt bs 


convenient to discuss them separately in the first instance. 


wealeotaqne irc | Preitd 


eo? Gee awe ae ee ae Oe ee 


It is common practice to average paleomagnetic data in 
a hierarchical fashion (1rVvVingd,~ 19653) Metlhinny, 1973) using 
the procedure due originally to Fisher (1953). . Specimens 
from.a single site (horizon) are considered to have acquired 
their primary nagnetization Simultaneously. Averaging 
specimen directions to obtain site means reduces the effect 
of experimental error associated with collection, 


orientation and measurement. The site means are then 
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combined to obtain the mean paleomagnetic direction for the 
entire section. This procedure provides estimates of the 
magnitude of different types of errors. The within-site 
scatter reflects experimental error and _ the recording 
fidelity of the strata, whereas the between-site scatter 


obviously includes the effect of real temporal variations. 


' The site means for the 0, 10 and 20 mT demagnetization 
data are summarized in Table 4.1 and illustrated in Fig. 
4{.1la-c. The overall mean direction for the entire sequence 
at O mT is D=006.39, I=+66.49 (N=94, R=93.3895, k=152, 
alpha95=1.29); at 10 mf it is D=005.79°, I=+65.49 (N=94, 
R=93.3099, k=135, alpha95=1.3°) and at 20 mT it is D=005.99, 
I=+64.99 (N=94, R=93.2860, k=130, alpha95=1.3°). All these 
are Significantly different from the present geomagnetic 
vector at the sampling locality (D=022.3°9, I=+72.6°) (see 
Fig. 3.1a). Since there is very little difference between 
the results at 10 and 20 mT, the site means at 10 mT will 
henceforth be used. On the whole, within-site grouping is 
good, with values of Fisher's precision parameter (k) 
greater than 100 in 85% of cases (78 out of 94) (minimum=20, 
maximum=5,242, arithmetic mean = 342). There is no 
correlation between k and the lateral spread of the 
specimens at each site (Fig. 4.2). This implies that the 
scatter is not due to non-horizontality and/or irregularity 


of bedding, or to inadequate vertical positioning of 
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Table 4.1 Site means at 0,10 and 20 ot. 
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Specinens; Le is probably due to minor lithological 


variations. 


The standard two-tier analysis (Watson & Irving, 1957) 
yields within-site (ky) and between-site (kg) precisions of 
140 and 178 respectively. As expected, the heey onieeres 
precision increases when freed of the effect of within-site 
scatter. Completion of the two-tier analysis leads to 
ko= 12,716 and alpha95=1.2°, but the slight improvement is of 


no real significance. 


The pole obtained from the mean direction lies at 
3. 4°ER, S5.4°N (dp=1.7°, dm=1.2°), whidst .the mean of the 
virtual geomagnetic poles (VGP's) (Fig.4.3) lies at 0.0°E, 
85.6°N (N=94, R=92.7032, K=72, A95=1.7°). Thus the overall 
mean paleomagnetic pole is statistically distinct fron the 


rotation axis of the earth. 


In early paleomagnetic work this small (4.49) 
diitference would probably have been attributed to 
methodological inadequacies, but a great deal of interest is 
currently being directed towards*® these ‘second order * 
features, and further scrutiny 4s therefore Warranted. fhe 
pole obtained is '‘'right-handed' and "Wal s10ecr, in 
agreement with global analyses of Wilson (1970, 1971, 1972). 
Merrill & McElhinny (1977) have updated these studies and 


provide a comprehensive compilation of 266 high-quality 
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Fig. 4.3 Virtual geomagnetic poles (VGP's ) for all 94 horizons 


with their. mean and its circle of 95% confidence. The 
dipole axis for the 1965 International Geomagnetic 
Reference Field and the pole corresponding to the 
present field vector at the sampling site are indicated 
by IGRF and PF respectively. The longitude of the 
sampling site is shown by an arrow at the perimeter. 
( Equal area projection ) 
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published results. Of these 101 are of particular interest 
to the present discussion, being of normal polarity and 
Quaternary age (0- 2 msy.). For discussion of the 'right- 
handed*® and 'far-sided' effects Wilson (1972) introduced the 
convenient ‘common site longitude* (CSL) concept in which 
each sampling site is imagined to be at zero longitude. 
Merrill and McElhinny (1977) obtained a mean CSL pole 
position of 87.72°N, 166.14°F (N=101, R=100.03221, K=103.3, 
A95=1.39). By this same convention, the Riggins Road pole 
lies at 85.6°N, 118.9°F with A95=1.7° (Fig. 4.4). These two 
poles differ by only 3.39; assuming a Fisherian distribution 
with K=103.3, the probability of observing such an 
individual divergence is quite high (84%). It thus appears 
that the Riggins Road data typify the time-averaged 
paleomagnetic field, which is dominated by the axial dipole, 
but has detectable higher order terms. This, in turn, 
implies that under suitable sedimentological conditions a 
tine spar of 9,000 years is sufficient to yieid an 
acceptable estimate of the local paleofield. However, poles 
deduced from such data using the usual dipole mapping 
function will obviously not yield a precise estimate of the 
spin axis (geographic pole), because of the persistent non- 
dipole terms. Merrill & McElhinny propose a slightly 


modified procedure for calculating poles which allows for 


these effects. 
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and for the Riggins Road data with 


respective circles of 95% confidence ( see text ). 


(engue.! yy .ared projection ) 
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Cox (1975) has summarized a large body of excellent 
paleomagnetic data from basaltic flows exposed on the 
Hawaiian Islands. He concludes that the time-averaged 
paleofield diverges from that of a geocentric axial dipole 
field by 7.7°+29, the discrepancy being entirely in 
inclination. The Smallest long-term ‘anomalous’ field 
Capable of explaining this offset. is) 5,400 nT ( AX=3,300) nT, 
AZ=-4,400 nT), which is definitely appreciable compared to 
typical IGRF non-dipole values. As possible explanations Cox 
rejects Wilson's (1971) offset dipole model, and also the 
‘standing! non-dipole field model of Yukutake & MTachinaka 
(1969). His favoured explanation involves biassed drifting 
fields and is discussed more fully in section 4.3 in 
connection with the present data. For the moment we note 
that the Riggins Road pole lies 4.4° from the spin axis, 
corresponding to a directional divergence at the site of 
3.19 The minimun long term anomalous field capable of 
explaining this offset is 2,800 nT, (assuming an axial 
dipole! ftela “of / 52,000 (nT ate 50.5°N) | with componencs 
Ax=1,700) ut, JAY=2, 100 nT, and’AZ-77 00 nt. 2tlseseenithat, 
unlike the Hawaiian case, the eastward (i.e. AY) component 


is very appreciable. 


Merrill & McElhinny (1977) feel that the permanent 
guadrupole is now sufficiently well established (at least 


for the last 5 my.) to warrant a modification such that the 
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usual ‘dipole equation' tanI = 2cotP is replaced by : 

tanI = [4cosP + 3R(cos2P-1) J/[2sinP + 3Rsin2P] 
Where I 1s inclination, P is colatitude and R=g0/g9 : g9 and 
g° being the geocentric axial quadrupole and dipole Gaussian 
coefficients respectively. This modification refers to the 
far-sided effect, the proposed field being entirely zonal. 
They obtain R=+0.050 for normal polarity data spanning the 
last 5 meye; this compares favourably “with the 1975 IGRF 
value of +0.063 . Taking R=+0.050 the corrected Riggins Road 
pote shespatyies. 190, 86.3°N. (ig. o4. 4), stil l-3.7° from the 
Spin axis and now requiring a ‘permanent anomalous field of 
2,200 nT (AX=400 nf, AY=2,100 nf, AZ=-200 nT). The eastward 
component now dominates of course, since the meridional 
components have been essentially removed by the quadrupole 


correction. 


Merrill & McElhinny also analyse the right-handed 
effect first noted by Wilson (1972). They conclude that 
Wilson was misled by unevenly distributed data, there being 
a paucity of data from the Pacific region. They propose both 
right-handed and left-handed anomalies constituting a 
symmetrical sinusoidal variation of peak amplitude +2° (Fig. 
4.5). At ‘the Longitude and latitude of Riggins Road they 
predict a left-handed divergence of 2.8° whereas we find a 
right-handed value of AD=5.7°. Obviously a single datum such 


as this. cannot rule out their nodel, but our data, which 
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longitude ~ sectors. 
AD, = tan“ ( tanAD cos) where AD is the 
declination anomaly ( departure of D from zero ) 
ano, pA else the, Jatiindes ihnusi. AD. is the 
equator-normalized value of AD. Simple sinusoidal 
curves that best fit thee. orn): data lie 
within . the hatched region. Ears i ‘indicate the 
standard error ( after 


Merrill & McElhinny, 1977 ). 
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seem to represent a high quality recording, clearly conflict 


with it. 


These differences are small and further work is 
necessary to fully establish their validity. At the present 
time the observations and their implications are simply 
presented at face value. However, the detailed analysis and 
discussion of the temporal sequence of field vectors given 
in the following section lends credence to the above 


inferences. 


4.3 Declination and Inclination Oscillations 


The geomagnetic record (Fig. 4.1a-c) reveals no 
evidence of any large angular shifts which could be 
correlated with the geomagnetic excursions reported fron 
Mono Lake, California (25,000-24,000 yrs.B.P.;Denham, 1974; 
Liddicoat & Coe, 1975) or that from Lake Mungo, Australia 
(30,780-28,140 yrs.B.Pe; Barbetti & McElhinny, 1972) even 
though the sampling interval iS approximately 100 years. 
This implies that such excursions are not global events, but 
result from the spatially restricted effects of specific 
non-dipole sources (Harrison & Ramirez, 1975; Coe, 1977). No 
large excursions are recorded, but smooth declination and 
inclination oscillations are observed, These presumably 


result from ‘regular’ secular variation and thus provide an 
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opportunity to study temporal geomagnetic variations. 


Beran & Watson (1967) have devised a test for serial 
correlation of a sequence of unit vectors. If N successive 
points are represented by directions of magnetization 
Xj, geeooogXy and if there is a strong nearest neighbour 
correlation, then the angle between successive vectors will 
be much smalier than the angle between vectors with very 
different serial numbers. Thus the cosines of the angles 
between successive vectors are expected to be large. 
Therefore a large L, where Li( 1) =>) X eee Re, Signifies serial 
correlation. The significance level is given by: 

SPHhl L(t) EL) eleva) jee 
where E(L) is the theoretical mean of L(P)=))X;.X,; for a 
completely randomized sequence of vectors and Var(L)=E({L) - 
E(L2). For the 1% significance level the value of s is 
2.326. For our data this test yields a value of s=8.860 
indicating that the probability of sequential ordering ‘is 
very bucy greater than 99%. This implies that the apparent 
oscillations in declination and inclination magnetograms 
(Fig. 4.1b) represent real temporal fluctuations in the 
paleofield. Similar oscillations were reported many years 
ago by Johnson et al. (1948), but only in this decade has 
there been a revival of interest in the study of secular 
Variation obtained fron recent sediments (Opdyke et 
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The oscillations are very clearly revealed when the 
data of Fige 4.1b are smoothed by a moving average (Fig. 
4.6a,b)- Fluctuations up to about 40° in declination and 
about 209 in inclination are observed, and visual inspection 
suggests periodicities of a few thousand years. Spectral 
analysis results will be more fully discussed “in the 
following section; at this stage it is both convenient and 
worthwhile to scrutinize the actual pattern of paleofield 
vectors. In Fig. 4.7a the sequence of 94 site mean 
directions are illustrated. A coherent pattern is 
discernible, but smoothing by a moving average emphasizes 
the underlying pattern (Fig. 4.7b,c). A series of open loops 
emerge, reminiscent of historical observatory records. fhe 
amplitude of these loops is similar to _ the historical 
records for London and Paris (see Fig. 1.2), but the 
periodicities involved are considerably larger. Perhaps the 
most remarkable feature is the fact that the same ‘three- 
quarter loop' path is traversed repeatedly, first ina 
counter clockwise sense, then clockwise, and finally half of 
it is traversed counterclockwise again. We return to the 
implications of this pattern below, but first it is 
instructive to compare these results with Cox's suggestion 
of biassed drifting fields introduced in section 4.2. The 
inset in Fige 4.7a (from Cox, 1975) indicates how he 


envisages spatially non-isotropic disturbances leading to a 
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Fig. 4.6a,b Magnetograms of three-point and five-point 
moving averages . 
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Bie hee Ce) Azimuthal equidistant plot of the 94 site means. 
@ oldest horizon 
m youngest horizon 
A vectorial mean of 94 sites 
¥ present field at sampling site 


O.. axial dipole direction 


The inset diagram (after Cox, 1975) illustrates 
idealized secular variation patterns (see text) 


(b),(c) Three-point and five-point moving averages of the 
data’ shownnin a) % 
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small offset in the time averaged paleofield vector. The 
Riggins Road data conform rather well to this concept, but 
the bias is strongly eastward as well as towards shallower 
inclinations. The overall picture is therefore quite 
satisfying, in terms of the final effect, but the underlying 
cause is not at all clear. Further discussion of the 
possible cause, or causes, is postponed until after the 


spectral analysis results are described. 


Fourier analysis of the declination and inclination 
Ragnetograms has been carFried out. The results (Fig. 4. 8a,b) 
confirm that the main periodicities are a few thousand 
years, as suggested by visual inspection. The power spectra 
obtained compare favourably with those of Yukutake (1962) 
who analysed Quaternary data from Japan and New England. 
Fourier analysis has several drawbacks for this type of 
short noisy data. An alternative method of spectral analysis 
whereby information concerning the degree and sense of 
looping can be obtained in addition to the power at any 
given frequency has been put forward by Denham (1975). The 
method essentially involves treating both the declination 
and inclination simultaneously as a complex time series, and 
uses the maximum entropy method (MEM) developed by Burg 


(1967). The spectra produced by complex series are often 
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Fig. 4.8a,b Declination and inclination power _ spectra 


obtained by Fourier analysis . 


yt gene ehereah ne _— 4 < Php ie he ta 
: oe 


a ie | 

on 

ie? 

te 
i~ Ps 
eS 
=e 

ot 


- rw ; ; Th . Al 
f bil ey eee etd a 
; i / Ph id 2h A 


4 + 


$ y \ ‘ ys ; 2 .f 


% « 4 rl ned 


iM an 
reg! / site 
Rose 
> i) 


ne : 


a a 


87 


asymmetrical ise. a peak in the 'positive' half of the 
spectrum is not the same size as its corresponding nae in 
the ‘negative’ half of the frequency spectrum. At any given 
frequency, the greater the degree of asymmetry, the more 
circular will be the looping. The looping at any arbitrary 
frequency, f , is generally elliptical with semi-major and 
semi-minor axes given by : 

Pee jee eee (st) "2 | and ICP (£) ]472-[ P(-£) J2724 
respectively where P(f) is the power associated with a 
Spectral peak centred on f. The sense of the looping is 
geterminea iby , the “sign of P(e) =P: (-f), (positive for 
clockwise and negative LOL counterclockwise). The 
appropriate complex series is generated by first mapping the 
directional information onto an azimuthal equidistant polar 
projection centred on the vectorial mean of the data. The 
MEM algorithm employed was originally designed to analyse 
real data only (Anderson, 1974) and minor modifications were 
necessary for analysing the complex time series; these were 
accomplished by Oberg (1978). Oberg@carried out a seriess of 
numerical experiments to test the method. He used a variety 
of synthetic data sets with known frequency content and 
noise. He found that the MEM technique is able to recover 
the input signal even in the presence of noise almost twice 
its anplatude, It Was therefore decided” toprappiy athe 


technique to the present data, using programs kindly 
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provided by C.J.Oberg. 


The ability of the maximum entropy technique to provide 
a good spectral estimate depends on the proper choice of 
prediction error filter (PEF) length. Ulrych and Bishop 
(1975) suggest that the optimum filter length is that which 
minimizes the so-called ‘final prediction error! (FPE). Fig. 
4.9 illustrates the FPE versus the filter length for the 
present data set. The position of the first mininum is 
rather surprising since a filter length of 3 would be far 
too short to resolve any appreciable detail in the spectrum. 
By using a filter length of 8, which is the second FPE 
minimum, a rather poor resolution of the spectrum is 
obtained (Fig. 4.10a). In order to obtain better resolution, 
filter lengths of 14 and 30 were also used (Fig. 4.10b,c). 
In all cases periods of about 6,000 and 2,000 years emerge 
(Table 4.2). These periods are similar to those found in the 
reconnaissance study by Oberg (1978), but there are some 
important differences. In particular, Oberg found a 
clockwise circular looping with a 2,000 year period, whereas 
the more complete data now available indicate both clockwise 


and counterclockwise looping at periods near 2,000 years. 


An attempt has been made to assess the frequency 
content as a function of time by subdividing the present 


data set. Fig.4.7a-c offers guidance as to the most 
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Bie. 4.9 Plot of FPE versus prediction error filter 
length . Ulrych & Bishop (1975) have suggested 
that a cutoff of N/2 (N = numer of data 
points) in the filter length be imposed when 
using the criterion . 
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The units associated with the ordinate can be 
thought of as (degrees)/unit frequency since the 
amplitude of oscillation (in degrees) associated 
With a “piven peak is equal to the square root 
of the area under the peak. Power inthe. positive 
frequency of the spectrum is indicative of clockwise 
looping whereas counterclockwise power resides in the 
negative frequency half of the spectrum . 
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Table 4.2 Periods and sense of looping with different filter 


lengths 


filter length 


sense of loop 
counterclockwise 
elockwise 


clockwise 


clockwise 


counterclockwise 


counterclockwise 
clockwise 


counterclockwise 
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appropriate manner in which to do this. A subdivision into 
three sections is suggested, ise. horizons 1 (elevation 6.71 
m) to 41 (elevation 3.91 m), 41 to 65 (elevation 2.21 mn), 
and 65 to 94 (elevation 0 m). MEM analysis of these three 
sub-sets yields the results illustrated in Fig. 4.11a-c and 
Summarized in Table 4.3. These confirm expectations based on 
visual inspection of Fig. 4.7a-c. The earlier part of the 
sequence is dominated by periods around 5,000-6,000 years, 
but the upper Bait is dominated by shorter periods of 2,000- 
3,000 years. This kind of variation in spectral content 
renders the se. of spectral estimators quite hazardous 
unless very long data sequences capable of meaningful 
subdivision are available. This is not usually the case in 
paleomagnetic investigations. Nevertheless the existence of 
major periodicities at 5,000-6,000 years and 2,000-3,000 
years seem to be well founded. Furthermore, it is worthwhile 
to note that despite the closer sampling used in this study 
compared to Oberg's reconnaissance work, no strong high 
frequency spectral peaks emerge. Reference to historical 
data (see Fig. 1.2) would suggest periodicities of several 
centuries. Although small peaks do appear in the MEM power 
spectra at these frequencies, they are too small to be 
regarded as real. This suggests that either such signals 
were absent, or that the deposition and dewatering involved 


in the natural recoding process somehow acts as a high-cut 
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Fig. 4.lla-c MEM power spectra of the divided data set. 
The filter lengths are chosen according to 
the Ulrych & Bishop (1975) criterion. 
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data set 
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filter and suppresses such signals. This lack of high- 
frequency content is supported by the results of the finer- 


scale carving experiment illustrated in Fig. 2.6. 


4.5 Discussion 


Some simple preliminary modelling is attempted in order 
to investigate the possible underlying cause(s) of the 
secular variation pattern observed. The non-dipole has been 
observed to undergo westward drift during the last few 
decades (Bullard et ale, 1950). At latitude 50°N this drift 
would produce the secular variation pattern of Fig. 4.12. 
Some features are similar to those observed in our data but 
the correspondence is not compelling. Furthermore, the 
periodicities would almost certainly be too short, unless 
the ‘usual* westward drift rate of about 0.2°/yr was a 
factor of two or three slower in the past. Periods of about 
5,000 years are more likely associated with the main axial 
dipole (Table 1.1). The model of Opdyke et ale (1972) 
involving a 6,000 year dipole precession is attractive. At 
latitude 50°N such a model would produce elliptical secular 
variation patterns (Fig. 4.13) which are similar but not 
identical to those shown in Fig. 4.7a-c. A model whereby the 
pole is allowed to ‘nod*' back and forth along a given 
meridian yields the temporal variations shown in Figs. 


4. 14a-e. These simple models show some of the required 
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Fig. 4.12 


Secular variations pattern obtained at 
latitude 50°N by Setting the 1945 
non-dipole field drifts westward (data 
fron Bullard et al., 1950) . 
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Fig. 4.14 (a) Temporal variations of inclination when the 
pole is alloweds to ‘"nod® back and Forth 
along the meridian of the observer (declination 
will always be zero) . (Numbers denote maximum 
displacement from the ‘spin axis in degrees ) 
(b),(c) Temporal variations of declination and 
inclination when the pole is allowed to ‘nod! 
back and forth along a meridian at 45° to 
the meridian of the observer . 
(d),(e) Temporal variations of declination and 
inclination when the pole is allowed to ‘nod! 
back and forth along a meridian at 90° to 


the meridian of the observer. 
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features but none of them fit the observed secular variation 
pattern at all well. We are therefore compelled to consider 
more sophisticated models. A variety of such models has been 
investigated and an example of these is illustrated in Fig. 
4.15. This model involves dipole oscillations (period 8,250 
years; see Fig. 1.3), dipole nodding along the geridtan of 
the observer (at 50°N) with amplitude 10° (period 5,500 
acenes and non-dipole oscillations (period 2,350 years). The 
non-dipole is assumed to be biassed eastward (maximum +5,000 
nT, minimum -1,000 nT). Such a bias is reminiscent of a 
"standing' non-dipole field as suggested by Yukutake 6& 
Tachinaka (1969) (Fags!) 4216)le lasttehes; context vit is worth 
noting that the present non-dipole field at the sampling 
Site also points eastwards (see Fig. 1.1). The similarity 
between the observed smoothed secular variation pattern and 
the model above is satisfactory for the preliminary 
purposes, and indicates that minor modifications could 


cLeproduce the observations very closely. 


Finally etek piles perhaps worth noting that Weyes (1978) 
has recently proposed that the spin axis ‘nods' with a 
period of 5,600 years. He bases his conclusion on analysis 
of ancient shorelines and proposes a link with glacial 
periods; the geographic pole movement amounts to only a 
degree or so. Although the connection with our geomagnetic 


period at 5,000-6,000 years is obscure, the similarity of 
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Fig. 4.15 Azimuthal equidistant plot of secular variations 


pattern obtained from a model whereby the 
Observer is situated “aur @acitude 50°N and 
(i). the dipole oscillates with a period 


Tho = 83250 years. 


(ii) the dipole nods along the meridian of 
the observer with a period Ton a 56500 
years. (dhe noddingweis assumed to bes. a 
Simple harmonic motion centred at the 
rotation axis and amplitude of 10°. The 
displacement from the rotation axis along 
the. .@iven weridian is) “civen\ by : 

X = 10sin(2mrt/Tp Ne 
Thus Sas ae vase is given by: 
= (M/R3) (1 + 3cosep)i/e op, 
ey Me (dipole moment) = (8 + Heosgrt/Tho) x 
10°2 A.m@ 
R (Earth's radius) = 6.356 x10° m. 
p paleomagnetic colatitude 
t time (years) 

(iii) the non-dipole is biassed eastwards and 
oscillates with ar period. img, = 2,350 
years. ‘The’ non-dipole antensity “is given 
Dy. 

‘ Fp = 3, 000sin(2rt/Typo) 42,000 per: 
Dashed lines represent the observed secular 
variations pattern (smoothed by five-point 
moving average ) (see Fig. 4.7c) 
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periods is rather intriguing. 
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CHAPTER 5 


Summary and Conclusions 


Nagata a al. (1943, 1945, 1949) studied the 
magnetization of the upper 7 metres of the Narita Hed which 
consists chiefly of fine sands spanning from 70,000-42,000 
yrSs.BeP. Yukutake (1962) smoothed the declination and 
inclination variations observed. Skiles (1970) observed that 
the smoothed data clearly define three clockwise loops with 
periods of 1,800 years, but the dating information is rather 
uncertain. Following a comprehensive discussion, Skiles 
concludes chat this Wii tpless a mewescNa rd. drift. 2oL the 
perturbing magnetic sources; he thus infers that westward 
drift of the non-dipole field was in existence some 50,000 


years ago. 


Opdyke et al. (1972) studied the paleomagnetism of two 
Aegean deep-sea cores together spanning from 27, 000-0 
y0SeB.P. The cores have average magnetic inclinations of 
5594149 sod. and 4694139 sede which are similar to the 
inclination of the geocentric axial dipole model at the 
sampling locality (I=+55°). No declination information is 
available, but the inclination variations reveal the 


presence of a 6,000 year periodicity. They interpret this 
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pattern in terms of ‘dipole wobble' in which the axis of the 
geomagnetic dipole precesses uniformly around the spin axis 


along a counterclockwise circular path of radius 12°. 


Mackereth (1971) describes long-period declination 
oscillations in cores taken from post-glacial organic 
sediments deposited at the bottom of Lake Windermere, 
England. Subsequent analysis by Creer et al. (1972) 
indicates that the corresponding inclination variations were 
of irregular period and amplitude. Thompson (1975) after 
extending these studies to other JK and Swedish lakes argued 
that the main characteristics of declination and inclination 
Variations of Lake Windermere are recognizable elsewhere in 
north-western Europe throughout the Holocene. It thus 
appears that paleomagnetic studies of Late Quaternary 
lacustrine sediments may have some potential for application 
to regional geological correlation, as well as providing 


essential data for geomagnetic studies. 


Creer (1974) obtained a paleomagnetic record from a 
core at station 1474 during the Black Sea cruise of 
AtlantisII of the Woods Hole Oceanographic Institution. Five 
radiocarbon ages were determined from samples taken at 
intervals along the 11 metre long core, which spans from 
25,000-7,000 yrs. BP. The magnetic record exhibits 


inclination variations with clearly defined maxima and 
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minima, and a period of about 2,800 years. Similar 
inclination oscillations were measured in a section of 
sediments from the Upper Dry Cave at Jeita, Lebanon spanning 
approximately 16,000-0 yrs.B.P. Creer & Kopper (1976) have 
correlated these records for the overlapping time interval. 
Thus the inclination oscillation appears to have continued 
from 25,000-0 yrs.B.P. The Jeita declination record yields 
evidence of a 1,700 year periodicity. The geomagnetic record 
from sediments of Lake Michigan (core stations 895 and 937) 
and Lake Eire (core station 13194) (Creer et al., 1976) are 
superficially similar to those obtained at Lake Windermere 
with distinct declination oscillations, but there are no 
strong periodicities in inclination. However at Lake 
Michigan, the declination is observed to have oscillated 
back to 11,500 yrseB.P. with a period of about 2,000 years 
which is somewhat less than the 2,700 year period observed 


at Lake Windermere. 


It thus appears that regular declination oscillations 
are present at some sites whilst at other sites inclination 
oscillations are more prominent. In both cases periods of a 
few thousands years are characteristic. The absence of a 
steady inclination period related to the declination period 
will mean that the declination oscillations cannot be 
accounted for simply by supposing that the geomagnetic axis 


precessed about the spin axis. Furthermore, when the path of 
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the time sequence of paleomagnetic field vectors is plotted 
on the "D-I" plane, both clockwise and counterclockwise 
loops are observed (Creer et al., 1972). Thus the observed 
Variation may originate from westward and sometimes from 
eastward drifting sources (Skiles, 1970). A model whereby 
the geomagnetic sources drift either westward or eastward 
around the geographic axis would produce similar secular 
variation patterns for sites edad approximately the 
same latitude, with a time lag due to the longitudinal 
difference. The observations, however, indicate periods 
which remain constant over many oscillations at a given 
site, and different dominant periods at different sites. 
Thus Creer (1977) favours a model in which the non-dipole 
field originates from stationary oscillating or pulsating 
sources situated in the outer core. The corresponding 
perturbing fields therefore have a regional rather than a 


global influence. 


The present work is a follow up detailed study of the 
reconnaissance investigation by Oberg (1978). A 7 metres 
thick section spanning some 9,000 years which is part of a 
much thicker sequences laid down during the Olympia 
Interglaciation (43,000-20,000 yrs.B.P.) was sampled at 7 cm 
intervals. Magnetic measurements reveal the presence of a 
small viscous remanent magnetization which is removed by AF 


demagnetization in 10 mf, leaving a stable remanent 
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magnetization. IRM curves and thermal demagnetization 
indicate that the magnetic carrier is titanomagnetite, and 
this suggests that the paleomagnetic directions represent a 
DRM rather than a later CRM which would normally be ae peerad 
to reside in hematite. The overall mean direction for the 
entire sequence is D=005.79, T=+65.49 (N=94, R=93. 3099, 
k=135, alpha95=1.39), which is significantly different from 
the geomagnetic field vector at the Jaren Site (D=022. 3°, 
I=+72.6°). The mean of the virtual geomagnetic poles lies at 
O.0°E and 85.69N (N=94, R=92.7032, K=72, A95=1.79) being 
both *right-handed! and ‘far-sided*, and is statistically 
distinct from the spin axis of the Earth, The application of 
the small quadrupole correction suggested by Merrill & 
McElhinny (1977) moves the pole in such a way that it is no 
Ponger ~'far-sided", but At Sis eestrell right-handed. This 
conflicts with Merrill & McElhinny (1977) who suggest a 


left-handed (negative AD) anomaly at the sampling locality. 


The Beran-Watson test indicates that the serial 
correlation of the geomagnetic vectors is much greater than 
99%, and this implies that the oscillations observed are 
real temporal variations of the paleofield. Periods of a few 
thousand years are observed; but no periods of around 500 
years (as would be expected from historical records) are 
present even though the sampling density has a Nyquist cut- 


off of approximately 200 years. The absence of these high 
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frequency sigrals in the geomagnetic record is probably due 
to the deposition and dewatering involved in the natural 
recording process which somehow acts as a high-cut filter 


and suppresses such signals. 


The secular variation pattern agrees very well with 
Cox's (1975) model of biassed drifting fields, but the bias 
is strongly eastward as well as towards lower inclinations. 
This is consistent with Wilson's ‘right-handed' effect, but 


the origin of such a global feature remains obscure. 
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Appendix 1 
List of data used fn *“Ghapter-~ 2 
S site 
E elevation ( metres ) 
AF alternating field demagnetization ( mT ) 
D declination ( degrees ) 
L inclination ( degrees ) 


J magnetic intensity ( A.m? / kg ) 
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Appendix 8) 


Viscous remanent magnetization data 
t time in hours 
veru@ viscous remanent magnetization acquired ( A.m@ ) 
AF AF demagnetization ( nT ) 
vr“ viscous remanent magnetization, remaining 
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Appendix 4 


IRM build-up data 
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